Preprint: 1-6 (1998) 



FOSSIL IMPRINTS OF THE FIRST GENERATION SUPERNOVA 
EJECTA IN EXTREMELY METAL-DEFICIENT STARS 

Toshikazu Shigeyama 

Department of Astronomy, School of Science, University of Tokyo, Bunkyo-ku, Tokyo, 113-0033 Japan 
Research Center for the Early Universe, School of Science, University of Tokyo, Bunkyo-ku, Tokyo, 113-0033 Japan; 

shigeyamaQastron. s.u-tokyo. ac.jp 
Takuji Tsujimoto 
National Astronomical Observatory, Mitaka, Tokyo, 181-8588 Japan; 
tsuji@misty.mtk.nao.ac.jp 

(Received 8 July 1998; accepted 14 September 1998) 

Abstract 

Using results of nucleosynthesis calculations for theoretical core-collapse supernova models with various 
progenitor's masses, it is shown that abundance patterns of C, Mg, Si, Ca, and H seen in extremely metal- 
deficient stars with [Fe/H] ^ -2.5 follow those seen in the individual first generation supernova remnants 
(SNRs). This suggests that most of the stars with [Fe/H] ^ -2.5 were made from individual supernova 
(SN) events. To obtain the ratio of heavy elements to hydrogen, a formula is derived to estimate the mass 
of hydrogen swept up by a SNR when it occurs in the interstellar matter with the primordial abundances. 
We use [Mg/H] to indicate the metallicities instead of [Fe/H]. The metallicities [Mg/H] predicted from these 
SNRs range from ~ —4 to ~ —1.5 and the mass of Mg in a SN is well correlated with its progenitor's mass. 
Thus the observed [Mg/H] in an extremely metal deficient star has a correspondence to the progenitor's 
mass. A larger [Mg/H] corresponds to a larger progenitor's mass. Therefore, so called 'age-metallicity 
relation' does not hold for stars with [Fe/H] ^ -2.5. In contrast, the [Mg/Fe] ratios in the theoretical 
SNRs have a different trend from those in extremely metal-deficient stars. It is also shown that the 
observed trend of [Mg/Fe] can predict the Fe yield of each SN given the correspondence of [Mg/H] to 
the progenitor's mass. The Fe yields thus obtained are consistent with those derived from SN light curve 
analyses. This indicates that there is still a problem in modelling a core-collapse supernova at its beginning 
of explosion or mass cut. The abundance determination of O in extremely metal-deficient stars, that have 
not been done from observational analyses, are strongly desired to test the hypothesis that the elements in 
an extremely metal-deficient star come from a single SN event and to obtain reliable yields for SNe. 
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the lowest metallicity ([Fe/H]~ —4). 



1. INTRODUCTION 

Recent observations and the analyses imply that the 
abundance pattern of an extremely metal-deficient star 
with [Fe/H] (= log (Fe/H) - log (Fe/H) Q ) < -2.5 may 
retain information of a preceding single supernova (SN) 



event c|r at most a few SNe ( Mc William et al. 1995 ; Ryan 



Norris, fe Beers 1996[ ). A theoretical a ttempt from this 
point of view has been already made (Audouze & Silk 



1995) to explain the observed abundance patterns. They 



argued that a combination of yields from two SNe with 
different progenitor's masses at the main sequence (M ms ) 
is consistent with the abundance patterns of stars with 



Assuming that the formation of extremely metal- 
deficient stars is triggered by a single supernova remnant 
(SNR) and that the formed stars retain the abundance 
pattern of this SN, one could predict abundance pat- 
terns (including hydrogen) of these st ars from theoretical 
SN models available at present (e.g., Woosley fe Weaver 
19951 ; [Tsujimoto et al. 1995[ |Nomoto et al. 1997|) . In this 



letter, we will compare the abundance pattern thus ob- 
tained with observations to test whether the chemical en- 
richment by individual SNRs could explain the observed 
abundances on the surfaces of extremely metal-deficient 



stars. We will use models in Woosley fe Weaver (1995 



(referred to WW95) and those in Tsujimoto et al. (1995) 



2 



T. Shigeyama & T. Tsujimoto 



and |Nomoto et al. (1997)| (T95) 



The metallicity of a star has been usually indicated 
by [Fe/H]. On the other hand, the yields of Fe from SN 
models to date have not converged (Woosley & Weaver 
1995|f rn sujimoto et al. 1995) because of uncertainties in 



the explosion mechanism and fall back dynamics or mass 
cut between the forming neutron star (or black hole) and 
the ejected envelope (e.g., Thielemann, Hashimoto, & 
Nomot |o"l990 ). Thus the yields of lighter a-elements are 
more reliably calculated. We will use [Mg/H] instead of 
[Fe/H] to specify the metallicity because of the following 



1. Mg is less affected by the mass cut in SN models 
than Fe. 

2. Mg is not synthesized or broken by the SN shock. 

3. The mass of ejected Mg increases with increasing 
progenitor's mass. 

4. The abundance of Mg is observationally derived for 
many stars with [Fe/H] ^ —2.5. 

There is a disadvantage in this element, that is, two SN 
models (WW95 and T95) give somewhat different masses 
of Mg as a function of progenitor's mass (squares in Fig. 
1). The best element to be used in this respect would be 
O (pentagons in Fig. 1) if O abundances were available 
for many stars with [Fe/H] ^ —2.5 (see McWilliam et 
al. 1995| ; [Ryan et al. 1996| ). 



. 8 



Mg 



M ms (M Q ) 



Fig. 1. The masses of Mg (circles) and O (pentagons) 
ejected from SN models are plotted as a. function 
of the progenitor's mass ( Woosley fe Weaver 1995| 
(open); Tsujimoto et al. 1995| (filled)). 



There is a clear trend in the observed [C /Mg] — [Mg/H] 
plot (crosses in Fig. 2) for [Mg/H]< —2 that cannot 
be reconciled with the conventional one-zone models of 
the Galactic chemical evolution that assume a complete 
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Fig. 2. The top panel: the crosses are the ob- 
served [C/Mg] for stars plotted against the [Mg/H] 
( [McWilliam et al. 1995j ). The open and filled cir- 



cles show the same quantities in the first genera- 
tion SNRs calculated from theoretical SN models 
( Woosley fe Weaver 1995| (open circles); |Tsujimotq 
et al. 1995 (filled circles)). The arrow indicates the 



change in the abundance pattern of t he model to 
reproduce th e SN 1987A observations ( [Thielemann 
et al. 1990 and references therein). 'The middle 



panel: Same as the top panel but for [Ca/Mg]. 
The bottom panel: Same as the top panel but for 
[Si/Mg], See the text for arrows. 



mixing of elements inside the Galaxy. The difference of 
the observed [C/Mg] in the range of -3.5 < [Mg/H] < -2 
gives a mean gradient A [C/Mg] /A [Mg/H] ~ -0.7. On 
the other hand, the one-zone Galactic chemical evolution 
model using T95 would predict the evolutionary change 
in [C/Mg] starting with the value of ~ —1.2 given by the 
most massive star M ms ~ 50Af Q and converging to the 
average value of ~ —0.9 integrated for stars in the mass 
range of 10 < M ms /M Q < 50 over the Salpeter initial 
mass function (IMF) with a slope of x = —1.35. It results 
in the corresponding gradient of A[C/Mg]/A[Mg/H] ~ 
+0.2 that has the opposite sign and a smaller absolute 



1 However, 



Audouzc & Silk (1995) 



have deduced that the chem- 
ical inhomogeneity is maximal in the range — 4 ^ [Fc/H] ^ — 2.5 
(-3.5 < [Mg/H] < - 2). 
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value than in observations. This argument also holds for 
WW95. If we introduce a metallicity-dependent IMF, the 
observed gradient requires a drastic change of the slope 

of the IMF from x 10 at [Fe/H] 4 to x ~ +0.5 at 

[Fe/H]~ —2.5. The index of x ~ —10 is almost equivalent 
to assuming that most of the stars with [Fe/H]~ —4 are 
formed from SNRs with M ma = 10 M©. It will be shown 
in the following section that this is indeed the case if we 
assume that most of the stars with [Fe/H] £ — 2.5 are 
formed from individual SN events. 

Here we will also predict the Fe yields as a function of 
M ms using the observed [Mg/Fe]-[Mg/H] trend combined 
with the [Mg/H]— M ms relation in theoretical SN models. 
Then these Fe yields are compared with those derived 
from SN light curve analyses. To calculate the mass of 
hydrogen swept up by a SNR, that is important to obtain 



[Mg/H ], the analytical expres sion for this mass in Ciofh 



Mckee, fe Bertschinger (1988) is modified to be applicable 
to SNRs that occur in a metal-free interstellar matter 
(ISM). 



2. LATE STAGE EVOLUTION OF THE 
FIRST GENERATION SUPERNOVA REM- 
NANTS 



In a spherically symmetric SNR (e.g., Cioffi et al. 1988), 
after the radiative cooling time scale becomes shorter 
than the dynamical time scale at the shock front, the 
shock front propagates considerably slower than before 
and makes a dense shell immediately behind it. Most of 
the mass inside the shock front resides in this shell. The 
ejecta that contain heavy elements occupy only a small 
central portion of the SNR even during this pressure 
driven snowplow (PDS) phase. In a realistic SNR, the sit- 
uation is different. First, the contact surface between the 
ejecta and the ISM is subject to the Rayleigh- Tayl or in- 



stability during the preceding Sedov- Taylor phase (Tay- 
lor 1950 ; Sedov 1959| ) . A large part of the ejecta pene- 



trates into the ISM. Thus it is expected that the ejecta 
already approach the dense shell at the beginning of the 
PDS phase and merge with it during this phase. Sec- 
ondly, an isothermal shock front is considered to be dy- 



namically unstable ( Vishniac 1983). Thus the dense shell 
behind the shock front is fragmented into a number of 
cloud cores that retain the abundance pattern of the SN. 
These cloud cores embedded in a high ambient pressure 
become seeds of stars of the next generation ( [Nakano 
1998). Some of the stars thus formed are currently ob- 
served as extremely metal-deficient stars. Accordingly, 
the average abundance pattern inside the SNR is assumed 
to represent that in all these stars. 

obtained the time t„ 



Cioffi ct al. (1 



Emerge 



when a 



SNR loses its identity. This time is a function of the ra- 
tio of the shock velocity to the sound speed Cs (or the 



velocity dispersion) of the ISM at the beginning of PDS 
phase (t = £pds) and ipDS itself. For a SNR in a homo- 
geneous ISM with the density of n\, ipDS is determined 
by equation 



1 de 
n\ dt 



= -ax A(T), 



(1) 



where e and A(T) are the internal energy and the cooling 
function at the temperature T, respectively. All the vari- 
ables are evaluated at the shock front using Sedov- Taylor 
solutions for the point explosion with the energy Eq and 
the number density n\. The constant a should be equal 
to 1.85 to reproduc e ipDS for the numerical simulation in 
Cioffi ct al. (1988) when A(T) for the solar abundances is 
used. Then h. v -[i m ord\a\{T) is constructed for the gas com- 
posed of only hydrogen and helium with their mass ratio 
Ah : Ahc = 0.75 : 0.25 under the collisional ionization 
equilibrium. For n% ~ 10~ 2 cm -3 , equation ([l]) gives T 
greater than 10 5 K. Thus the ionization of hydrogen does 
not affect the dynamics. 

The mass of hydrogen M sw thus obtained with A(T) = 
AprimordialCO m equation (Q) is approximated by the for- 
mula 



M sw = 5.1 x 10 4 M o 



( Eo \ 
\l0 51 cr g/ l 



0.97 



-0.062 



Cs 



lOkms" 



-9/7 



(2) 



The mass M sw is insensitive to n\. The sound speed 
C s was assumed to be 10 km s" 1 (or T ~ 10 4 K). This 
mass depends on Eq and n\ in a different way from that 
of Cioffi ct al. (1988) due to the different cooling func- 



tion used here. The cooling function with the primordial 
abundances is approximately proportional to T~ 2 instead 
of j"-i/2 in|Cioffi ct al. (1988)|. 



3. ABUNDANCE PATTERNS 

Here we assume that all the stars with [Mg/H]< —2 are 
made from individual SNRs. Abundance patterns in in- 
dividual SNRs will be calculated based on yields from 
theoretical SN models and compared with t hose on the 
surface of m etal-deficient stars obtained by McWilliam 
et al. (1995) . For SN models in metal-free environment, 
we use models Z12A, Z13A, Z15A, Z22A, Z30B, Z35C, 
Z40C in WW95 to retain the monotonicity of the mass of 
Mg as a function of progenitor's mass (see filled squares 
in Fig. 1). Another group of SN models are taken from 
T95 in which the initial metallicity is solar. Parameters 
of both models are listed in Table |l|. 

We divide elements into two categories. One category 
is composed of the elements not (or at least less) affected 
by the mass cut in SN models. The examples are C, 
Mg, Si, and Ca. The other is composed of those affected 
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Table 1. Input parameters in SN models and the calculated [Mg/H] ratios in the corresponding SNRs (see text for details). 





WW95 






T95 




Name 


M ms (M ) 


E Q (10 51 erg) 


rA/Tn- /"Hi 

[Mg/ ti\ 




-Bo (.10 erg) 


TA/Tn- /"Hi 

[Mg/ i±\ 


Z12A 


12 


1.28 


-4.1 


13 


1 


-3.6 


Z13A 


13 


1.29 


-3.6 


15 


1 


-3.1 


Z15A 


15 


1.27 


-3.2 


18 


1 


-2.9 


Z22A 


22 


1.26 


-2.8 


20 


1 


-2.4 


Z25B 


25 


1.83 


-2.9 


25 


1 


-2.3 


Z30B 


30 


2.06 


-2.5 


40 


1 


-1.9 


Z35C 


35 


2.49 


-2.6 


70 


1 


-1.6 


Z40C 


40 


3.01 


-2.6 









by the mass cut and includes Cr, Mn, Fe, Co, Ni, for 
example. This division is necessary because the mass cut 
in SN models to date is totally artificial and uncertain. 

Elements not influenced by mass cut — The yields of el- 
ements not influenced by the mass cut are expected to be 
more reliably calculated in SN models than those influ- 
enced by the mass cut. Therefore observed abundance ra- 
tios of [C/Mg], [Ca/Mg], [Si/Mg] with respect to [Mg/H] 
(Fig. 2) can tell whether extremely metal-deficient stars 
were formed from individual SNe. These abundance ra- 
tios derived from SN models (open circles:WW95; filled 
circles: T95) are also plotted in the same figure, in which 
[Mg/H] is calculated from the ratio of the mass of Mg in a 
SNR to that of hydrogen swept up by the same SNR (see 
equation (2)). The relations between [Mg/H] and M ms 
are shown in Table 1 for both SN models. The values 
of [Mg/H] for M ms ~ 10M Q (the lower mass limit of the 
core-collapse SN progenitor) from both models are not so 
far from the observed lowest value of [Mg/H]^ —3.7 and 
the metallicity [Mg/H] increases with increasing progen- 
itor mass M ms . This infers that there are more stars at 
lower metallicities according to an IMF that decreases to- 
ward high masses if stars with lower masses, say 10 M Q , 
explode in the regions that have not been polluted by 
other SNRs before. This condition is represented as 
™*Vsnr Jiom® <t>{m)dm < 1 star. Here n*, <p, and Vsnr 
denote the number density of the first generation stars, 
the IMF normalized to unity between the lower and up- 
per mass limits, and the maximum volume occupied by a 
single SNR, respectively. This condition results in a star 
formation efficiency for the first generation stars smaller 
than 2%, if the Salpeter IMF is used. 

A trend seen in the observed [C/Mg] (crosses in the top 
panel of Fig. 2) is well reproduced by both models. This 
means that stars with [Mg/H]~ —3.5 were made from 
individual low mass (M ms ~ 10 M©) SNe and that stars 
with [Mg/H]^ —2 were from individual high mass SNe 
(M ms ~ 50 Mq). For [Si/Mg], both models are also con- 
sistent with observations, though the observed [Si/Mg] 



has no clear trend. The [Ca/Mg] ratios derived from 
T95 with M ms > 20M Q deviate from observations and 
gives lower values in the region [Mg/H]> —2.5. This in- 
dicates that the yields of Ca of T95 with M ms > 20M Q 
are too small. These abundance patterns together with 
the predicted [Mg/H] ratios in Table 1 suggest that stars 
with —3.5 ^ [Mg/H] ;$ — 2 are made from individual SN 
events. 

Elements influenced by mass cut — The upper panel of 
Fig. 3 shows that the [Mg/Fe] ratios derived from both 
models do not follow the observations. This reconfirms 
that none of the two models correctly describes the dy- 
namics and/or nucleosynthesis near the surface of the Fc 
core in progenitor stars during explosion. In both models, 
less massive stars have too large Fe yields. More massive 
stars in T95 have too small Fe yields but the opposite is 
the case for WW95. The lower panel of Fig. 3 shows that 
the [Ca/Fe] ratios from both models are smaller than the 
observations. This can be ascribed to too large Fe yields 
from less massive stars of theoretical models again. An 
apparent fit of the observational data to the model for 
[Mg/H]> —2.5 is a consequence of two errors in yields 
of both Ca (see section ||) and Fe for more massive stars 
(see the next section). 

Abundances of some heavy elements in the ejecta of 
the 20 Mq supernova were deduced from the SN 1987A 
observations (Thielemann et al. 1990 and references 
therein). As the mass of Mg ejected from SN 1987A 
was not estimated from observations, we derive the mass 
of Mg by reducing the value of T95 to match the ob- 
served [Mg/Fe] ratio in extremely metal-deficient stars. 
The corresponding change of [Mg/Fe] is shown by the ar- 
row in the upper panel of Fig. 3. If we use the mass of Mg 
thus derived and the observed masses of C, Si, Fe, and Ca 
(the upper limit) , the abundance ratios of the model with 
M ms = 20 Mq would move following the arrows in Fig. 
2 and in the lower panel of Fig. 3. Thus the abundance 
pattern of the best observed SN coincides with that in 
extremely metal-deficient stars with [Mg/H]~ —2.6 and 
the mass of Mg in the 20 Mq model of T95 might be 
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but hardly change the average Fe yield integrated over 
the Salpeter IMF. Open circles are those from WW95 and 
give a similar amount of Fe for each progenitor's mass. 
Points with error bars also shown in the same figure are 
the Fe masses obtained from SN light curve analyses. 
Good fits of the Fe masses from the light curve analyses 
to those inferred from the abundance patterns demon- 
strate that this procedure to obtain the Fe yields works, 
if the Fe yields are not affected by the initial metallici- 
ties. A discrepancy of Fe masses at M ms ~ 10 M© might 
be due to the small number of fairly scattered observed 
points of [Mg/Fe] at [Mg/H]~ —3.5 in the upper panel 
of Fig. 3. The Fe yields thus obtained predict that the 
lowest [Fe/H] should be ~ —3.7 for both models. 



1 1 
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' 1 

98bw 



97ef 



1941 87A 
I 

93J | 



[Mg/H] 



Fig. 3. The upper panel: the crosses are the ob- 
served [Mg/Fe] for stars p lotted against the [Mg/H] 
(|McWiiliam et al. 1995| ). The solid curve shows 



the x- s< 3 uarc fi t of the observed points with the 
errors given in McWilliam et al. (1995) to a cu- 
bic polynomial formula . The open and filled cir- 
cles show the same quantities in the first genera- 
tion SNRs calculated from theoretical SN models 
("Woosley fe Weaver 1995 (open circles); Tsujimotc 



et al. 1995 (filled circles)). The lower panel: Same 
as the top panel but for [Ca/Fe], See the text for 
arrows. 



overestimated. 



4. IRON YIELDS DEDUCED FROM EX- 
TREMELY METAL DEFICIENT STAR OB- 
SERVATIONS 

We will reverse the argument in the preceding sections. 
Suppose that the trends in the observed abundance ra- 
tios, say [Mg/Fe], tell yields of SNe as a function of 
[Mg/H]. The x 2 fit °f the observational data to a cu- 
bic polynomial formula gives the solid curve in the upper 
panel of Fig. 3. This curve and [Mg/H] of each model 
predict the mass of Fe (Mp e ) from a SN as a function 
of M nls . Filled circles in Fig. 4 show the predicted Fe 
masses using T95. These Fe masses are completely dif- 
ferent from those given in the original models (triangles) 



M ms (M ) 



Fig. 4 The mas ses of Fe in SN models (jWoosley & 



Weaver 1995| (open); [Tsujimoto et al. 1995j (filled)) 
are shown by triangles. Circles show the masses of 
Fe in SN models with 1 a error bars of the x 2 fit m 
the upper panel of Fig. 3. These masses are deter- 
mined by the [Mg/Fe] — [Mg/H] relations indicated 
with the solid curve in the upper panel of Fig. 3 
and the [Mg/H]— M ms relations. Points with er- 
ror bars are the masses of 56 Ni, that eventually 
decays to 56 Fe derived from SN light curve p.na.l. 
yses (S N 1987A: Shigevama fe Nomoto 
1 9Q3.T 



. Shigeyama et al. 1 aa^ SN 1 HQ4T- ^ 

ct al. 19 94t SN 1997cf: jlwamoto et al. 19984 SN 



1998bw: 



-QQcf SN 
bgamotd 



Iwamoto et al. 



1998b) ). 



5. CONCLUSIONS 

We have shown that the abundance patterns of C, Mg, Si, 
Ca, and H theoretically predicted in the first generation 
SNRs are in good agreement with those observed in ex- 
tremely metal-deficient stars with —4 <[Fc/H]< —2.5. 
All these elements are thought to be less affected by 
the mass cut in SN modelling than heavier elements like 
Fe. This agreement implies that each extremely metal- 
deficient star was formed from a single SN event. In con- 
trast, both theoretical SN models predict different trends 
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in the [Mg/Fe] — [Mg/H] plot from the observed one. Con- 
versely, the mass of Fe ejected by each SN as a function of 
M ms can be derived from the observed [Mg/Fe] — [Mg/H] 
trend combined with the [Mg/H]— M ms relation in the- 
oretical SN models. This Mp e as a function of M ms is 
consistent with that derived from SN light curve analy- 
ses. Following the same procedure presented in this letter 
to obtain the Fe yields, we can derive the yields for other 
element s such as Ti, Cr, Mn, Co, Ni and r-process ele- 
ments ( Tsujimoto fe Shigeyama 1998| ), all of which arc 
very uncertain in SN models. If O abundances are de- 
duced in extremely metal-deficient stars, we will be able 
to obtain more reliable yields for SNe. In the forthcom- 
ing paper we will also discuss the connection between 
the chemical enrichment by individual SNe and the sub- 
sequent chemical evolution. 
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